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Abstract 

Born-Infeld black holes in the Scalar- Tensor Theories of Gravity, in the case 
of massless scalar field, have been recently obtained [H [21 [3]. The aim of the 
current paper is to study the effect from the inclusion of a potential for the 
scalar field in the theory, through a combination of analytical techniques and 
numerical methods. The black holes coupled to a massive scalar field have richer 
causal structure in comparison to the massless scalar field case. In the latter case, 
the black holes may have a second, inner horizon. The presence of potential for 
the scalar field allows the existence of extremal black holes for certain values 
of the mass of the scalar field and the magnetic (electric) charge of the black 
hole. The linear stability against spherically symmetric perturbations is studied. 
Arguments in favor of the general stability of the solutions coming from the 
application of the "turning point" method are also presented. 



1 Introduction 

The study of black holes coupled to non-linear electrodynamics is natural since in the 
conditions of strong fields and strong sources, such as black holes, quantum corrections 
to the lagrangian of electrodynamics should be taken into account. Nonlinear electro- 
dynamics was considered for the first time by Born and Infeld in 1934 in their attempt 
to obtain a finite energy density model for the electron The interest in nonlinear 
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electrodynamics has been recently revived since such types of lagrangians appear in 
the low-energy limit of open strings and D-branes [S] [H]- Nonlinear electrodynamics 
models coupled to gravity have been discussed in different aspects (see, for example, 
[9] [26] and references therein). 

The properties of charged black holes coupled to Born-Infeld nonlinear electrody- 
namics in General Relativity have been examined in [T3J [271 123 1.29] and in the case 
of non-zero cosmological constant in [30] . The role of the derivative corrections to the 
properties of the Einstein-Born- Infeld black holes have been studied in [31]. A natural 
step in the study of charged black holes is to add a scalar field in the theory. Born- 
Infeld-dilaton (asymptotically non-flat) black holes have been reported in [32] [39j . 
The case of charged black holes with massive dilaton coupled to abelian gauge field 
with linear lagrangian have been considered in [ID] - [12]. Einstein-Born- Infeld black 
holes with massive dilaton have also been investigated [^3] . 

In the frame of Scalar- Tensor Theories (STT) solutions describing charged black 
holes with non-trivial massless scalar field in non-linear electrodynamics have been 
recently obtained [H 121 E]. The aim of the current paper is to extend [TJ [3] by adding 
a potential for the scalar field that admits the presence of asymptotically flat black 
holes. 

The solutions with massless scalar field presented in [U [3] have a much simpler 
causal structure than the corresponding solutions in the Einstein-Born-Infeld (EBI) 
gravity. The latter have a single horizon, thus their structure resembles that of the 
Schwarzschild solution. Adding a potential for the scalar field makes the causal struc- 
ture richer and more difficult to study. In that case solutions with internal horizons 
and with degenerate event horizons may appear. 

A prior intuition about the properties of the studied solutions could be obtained 
from other similar systems that have already been studied. An example of such systems 
are the charged black holes with massive dilaton coupled to linear electrodynamics. Not 
all techniques used in that system are applicable to our case, though. The major dif- 
ference between the dilaton and the gravitational scalar in the STT is its coupling to 
mattei0, and in particular, to the electromagnetic field. In STT due to the specific cou- 
pling between the scalar and electromagnetic fields the source term in the field equation 
for the scalar field is proportional to the trace of the energy-momentum tensor of the 
electromagnetic field and not proportional to the lagrangian of the electromagnetic 
field as it is in the dilaton gravity. 

The charged black holes with massive dilaton have been studied analytically in [30] 
and numerical solutions have been obtained in [H]. In the first paper the conclusions 
about the causal structure are made. In that case black holes may have three, two 
or one horizons depending on the values of the parameters and the specific choice for 
the potential of the dilaton. The various types of extremal solutions were also studied. 
The authors found that two-fold and three-fold degenerate horizons may exist. The 
behavior of the fields near the central singularity was also obtained in that paper. In 
[H] approximate solution in the two limiting cases - large black holes (with radius 
of the event horizon much bigger than the Compton length of the dilaton) and small 

lr The coupling between the electromagnetic field and scalar field in dilaton gravity and in STT is 
different. The difference between the two types of actions can be found in |44j . for example. 
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black holes (with radius of the event horizon much smaller than the Compton length 
of the dilaton) have been obtained. Numerical solutions covering both the exterior and 
interior regions of the black holes are also presented. 

In the current paper our aim is to study the black holes solutions in a class of 
STT with massive scalar field coupled to Born-Infeld nonlinear electrodynamics. In 
the studied region of the parameter space we exclude the possibility of the existence 
of a third horizon. We also study the possibility for existence of extremal black holes 
in different regions of the parameter space and obtain numerical solutions for them in 
the allowed region. 

In the present paper we have commented also on the linear stability of the obtained 
black-hole solutions against radial perturbations. Here we follow the same scheme as 
for the black holes with massless scalar field [3]. Additional arguments in favor of the 
general stability of the black holes can be given also on the bases of the "turning point" 
method that we also discuss. 



2 Formulation of the problem 

The action of the studied STT is originally formulated in the so-called Jordan frame in 
which the scalar field is coupled to the scalar curvature and in which there is no direct 
coupling between the scalar field and the sources of gravity (in our case the source of 
gravity is the electromagnetic field). For mathematical convenience, however, as usual 
for the STT, we study the solutions in the conformally related Einstein frame (for more 
details we refer the reader to [IHl HHJ [1]). In the Einstein frame the action takes the 
following form 

S = y^q- J d *xV=9 [K ~ Wd^dvV - iV(<p)} 

+S m [^ m ;A 2 (ip)g^} (1) 

where TZ is the Ricci scalar curvature with respect to the Einstein metric g^ u and 
is the bare gravitational constant. 

The action of the nonlinear electrodynamics is 



l — J d*x^A 4 ( V )L(X,Y) (2) 



where 



X = ^^F tMU g^F a „ Y=^MF, u (*Fr (3) 



and "*" stands for the Hodge dual with respect to the Einstein frame metric g^ h 
The action ([T]) with §2§ yields the following field equations 
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U„ v = 2d^d^ + 2V{<p) gia , - 2d x L(X, Y) - ^g^F a(] F a ^ 

-2A\ V ) [L(X, Y) - Yd Y L(X, Y)] g, u , 

V„ [d x L(X, Y)F>» + d Y L(X, F)(*FT] = 0, (4) 

V„V>V = _ 4a( V )A\ V ) [L(X, Y) - Xd x L(X, Y) - Yd Y L(X, Y)] , 

where a(<p) = dln ^ Lp ' > ■ 

In what follows we consider the truncated^] Born-Infeld electrodynamics described 
by the Lagrangian 

L BI (X) = 2b I 1 - J I - I (5) 




The type of STT is determined by the specific choice of the function A((p) (re- 
spectively ot(ip)) and V(<p). Here we will consider potentials V(<p) which satisfy the 
conditions 

dV(<f) , dVUp = 0) 
<p , > and ^7 = 0. 6 

dip dip 

In particular, for numerical calculations, we will take the potential of the scalar field 
in the form 

V(<p) = (7) 

where m* is the mass of the scalar field and has also the meaning of inverse Compton 
wavelength of the scalar field in the units we work. Constraints for the values of m* 
can be found in [4"T] . 

In the present paper, we will be searching for solutions with regular scalar field <p> 
on the event horizon. We will also require that < A(ip) < oo for r > r#, where r# 
is the radius of the horizon in order to ensure the regularity of the transition between 
the Einstein and the Jordan conformal frames. Yet, we will consider only theories for 
which a(ip) has a fixed positive sign for all values of ip. The manner of investigation 
of solutions within theories with a(ip) < is similar. Theories in which the coupling 
function changes its sign are much more complicated (also from numerical point of view) 
since in them some interesting non-perturbative effects like "spontaneous scalarization" 
may appear, especially when a((p) ~ cp [2]. 

For our numerical solution we have considered theories for which the coupling func- 
tion has the form 

A(ip) = e«* (8) 

where a is a positive constant and in this theory az(ip) = const = a . We have 
studied the parametric space for fixed value of the coupling parameter a = 0.01. The 



2 Here we consider the pure magnetic case for which Y = 0. 
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conclusions for the qualitative behavior of the solutions, however, are valid for a much 
wider class of STT for which a(<p) > 0, a(<p = 0) ~ 10" 4 + 10" 2 and /% = 0) > 
—4.5 (such values are in agreement with the observations, for details see [IE]), where 
f3(<p) = da ^ . The exterior region solutions from the mentioned class of STT would 
not only be qualitatively the same but also quantitatively very close to the exterior 
region solutions in the Brans-Dicke theory since, as the numerical results show, ip is 
very small there and higher order (in (p) terms in the coupling function would have 
negligible contribution. 



3 Basic equations 
3.1 The reduced system 

In the present paper we will be searching for static, spherically symmetric, asymptoti- 
cally flat black holes. The metric of a static, spherically symmetric space-time can be 
written in the form 

ds 2 = g^dx^dx" = -f{r)e- 25{r) dt 2 + + r 2 (d6 2 + sin 2 Odcf) , (9) 

where 

/<r) = l-5=£. (10) 
r 

m(r) is the so-called local gravitational mass and should not be confused with the mass 
of the scalar field m*. We will study the magnetically charged black holes for which 
the electromagnetic field strength is given by 

F = P sin 9d9 A dcp (11) 

and the magnetic charge is denoted by P. The electrically charged solutions can be 
obtained through electric-magnetic duality rotations (we refer the reader to [TBI dH [TBI 
|3] for more details on these duality transformations) of the type 

{g^, <p 7 P, X, L(X)} <— > { 9llv , <p, *G^, Q, X, L(X)}. (12) 

The barred quantities are related to the dual system. Here 



X = - 



d x L(X) X, (14) 



and Q is the electric charge of the dual solution. The Born-Infeld lagrangian (jSJ) 
belongs to the class of lagrangians for which the system (J4J) is invariant under the 
electric- magnetic duality transformations (1121) . As a consequence of that invariance 
the solutions for the metric functions and the scalar field in the electrically charged 
case and the magnetically charged case coincide. 
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The field equations reduce to the following coupled system of ordinary differential 
equations 



/" - 2f5" - 3f'S' + 2f5' 2 + -f - -f6' = 

-4 {V(tp) + A\ V ) [2Xd x L(X) - L(X)}} , (15) 

/" - 2f5" - 3f'5' + 2f5' 2 + -f = 

-4 { V 2 + V(tp) + A\<p) [2Xd x L(X) - L(X)}} , (16) 
l-f-rf' + rfS' = 2r 2 [V(<p) - A(<p)*L(X)] , (17) 
d 
dr 



,dip\ = dV(<p) 
dr J dip 



- ( e"^ 2 /^ ) = + Ar 2 e- s a( V )A\ V ) [Xd x L(X) - L(X)} , (18) 



where X reduces to 



X = (19) 

It is a system of four equations for three unknown functions but, as it is well known, 
the self-consistency of the system is guaranteed by the Bianchi identity. For numerical 
treatment of the problem, however, the following form is more convenient 



dS ( dtp 

dr \ dr 



dm 
dr 



= r 



i/(-V 

2 J \dr J 



V(ip) - A(ipfL(X) 



dr 



dr 



dV 



— - 4a(¥>M» \ L ~ Xd x L{X)\ ~rf^ 



dip 



(20) 
(21) 

(22) 



In this form the last two equations are separated as an independent sub-system. 



3.2 Qualitative investigation 

Some general properties of the solutions can be derived through an analytical investiga- 
tion of the equations. The major difference resulting from the different type of coupling 
between the scalar field and the electromagnetic field is as follows. In both theories, 
STT and dilaton gravity, the right-hand side of some of the generalized Einstein equa- 
tions for the metric functions is proportional to the lagrangian of the electromagnetic 
field. The right-hand side of the equation for the scalar field in STT is proportional 
to the to the trace of the energy-momentum tensor of the electromagnetic field while 
in dilaton gravity it is proportional to the lagrangian of the electromagnetic field. The 
last fact allows in the case of dilaton gravity combinations between the field equations 
which are favorable for the analytical investigation of the properties of the solutions to 
be made. 

Analytical assessments and results, though not very complete, are always welcome 
since the entire space of parameters cannot be examined numerically thoroughly. The 
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examination of the properties of the solutions is additionally impeded by the presence 
of numerical instabilities for some values of the parameters. Such instabilities appear 
especially in the interior region and have been reported also in [41j . 

For our assessments we will use the fact that for the Born-Infeld lagrangian of 
nonlinear electrodynamics (3J the relations 



Xd x L(X) - L(X) > and 2Xd x L(X) - L(X) < 



(23) 



hold. 

For the analytical analysis we will also use the boundary conditions presented in the 
Section [5] where the numerical treatment of the problem is presented. More specifically, 
we will want the scalar field to be vanishing at the spacial infinity. 

It can be proved that the function (p(r) is non-positive outside the event horizon. 
Let us first exclude the possibility p(r) > on the event horizon. 

Proposition 1 p(r) cannot be positive on the event horizon. 

Proof: In order to prove this we will use equation ( 1T81) . Let us suppose that ip(r) is 
positive on the event horizon and has one or more roots in the exterior region. Then 
we multiply equation (118j) by (p(r) and integrate it in the interval r 6 [rjj, tq] where we 
denote the radius of the event horizon by r# and the leftmost zero of (p(r) outside the 
event horizon by r 



<P 



dr 



-<5„2 



r7 



dp 
dr 



dr 



rii 



ro 



dV(<p) 2 
<p — r 

dip 



+ Ar 2 e~ d p a(p)A\p) [Xd x L(X) - L{X)\ \ dr. (24) 



'7/ 



After integrating by parts we get 
pe r f 



dr 



ipe~ 5 r 2 f 



.dtp\ 



dr J 



ro r 



l-II 



I'll 



-5J2 



<hp\ 
dr J 



dr 



MS 



dr 



dV iv) 2 , A 2 - 

p> — r + 4r e 

dp> 



V a(ip)A 4 (ip) [Xd x L{X) - L(X)\ \ dr. (25) 



Taking into account (J6j) and (|23|) and the fact that according to the admission ip 
is positive in r e [^_f/> r o] we see that the sign of the right-hand side (RHS) of (1231) is 
positive. The left-hand side (LHS), however, is negative since /(r) > outside the 
event horizon. The contradiction we reach means that our admission is incorrect. If we 
admit that p(r) is positive on the event horizon and has no roots following the same 
procedure we will again reach a contradiction. In order to see this it is enough to let 
r — » oo in (j24"|) and to impose the boundary condition p(r) ^0 at infinity. □ 
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Proposition 2 Function (p(r) is negative and has no zeros in the exterior region. 

Proof: Let us admit that (f(r) has zeros. We know already that ip(r) is negative on 
the event horizon and also require it to satisfy the vanishing boundary conditions at 
infinity. Then one of the following two situations must be realized. Either the <p(r) 
has at least one positive maximum or it has at least one zero maximum. First, let 
us exclude the possibility for existence of positive maxima of if. Again, we can use 
( jMj) but this time let us integrate in the interval r e [Vq, r e ]. Quantity r e denotes the 
point in which ip(r) has a positive maximum, i.e., <f'{r e ) = 0, where (•)' denotes the 
derivative with respect to the radial coordinate r. Here r denotes the first zero of ip 
to the left of r e . After integrating by parts we get 



dr 



\dr J 

J d ^~ 1 ' 2 + 4rV< V ^)A 4 {cp) \Xd x L(X) - L(X)) | dr. (26) 



ro 



Since (p(r) is positive in r G [r , r e ] the sign of the RHS of (1261) is positive. The LHS, 
however, is negative. Again we reach a contradiction which means that the admission 
for presence of a positive maximum of <p is false. 

It remains only to exclude the case in which ip becomes zero at a maximum, which 
we will again denote as r e . At a zero maximum <p(r e ) = 0, (p'(r e ) = and f"{r e ) < 0, 
so from equation (1181) taking into account also (jSJ) we obtain 

fip"\ re = 4a(ip)A\y>) [Xd x L{X) - L(X)]\ re . (27) 

The LHS of the above expression is negative while the RHS is positive. This contradic- 
tion allow us to exclude the cases of zeros at maxima of (p. This completes the proof 
of the proposition. □ 

In our numerical calculations the following restriction in the behavior of / will be 
useful. 

Proposition 3 The metric function f(r) cannot have extrema with f(r) > 1. 
Let us consider equation (1171) . Using (1201) we can write it in the form 

l-f-rf- r 2 fp' 2 = 2r 2 [V{<p) - A(^) A L(X)} . (28) 
From here we obtain 

"1 



-r/'=(/-l) + 2r' 



2 f^ + V{fp) - A(<pfL(X) 



(29) 



When f{r) > 1 the RHS of (129]) is positive so the derivative of /(r) is strictly negative, 
f'(r) < 0. Hence, f(r) cannot have extrema with /(r) > l.D 

Once becoming greater than one, f(r) must continue to rise, no matter whether we 
go forward or backward in r. The inward integration can stop once the point f(r) = 1 
is surpassed. 
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Proposition 4 From the equations we can see that the function 5(r) is smooth non- 
increasing and has inflexion points where ip(r) has extrema. 



Proof: From equation (1201) we see that its first derivative is negative so the first part 
of the proposition is easily proven. Now let us differentiate equation fl20|) 

^ = _(^y_2r^.^. (30) 
dr 2 \dr J dr dr 2 

From (pOll and (!30il we see that in the extrema of <p(r) (the points in which ip'(r) = 0) 
5'(r) = S"(r) = 0, i.e., 5(r) has an inflexion point. □ 

3.3 Asymptotic solutions 

The presence of potential changes considerably the asymptotic behavior of the scalar 
field in comparison to the massless scalar field case. In the asymptotic region r — > oo 
we find the following behavior for the functions 

2 M P 2 P 4 

f(r) = 1 + — + — — + 0(r~ 14 ), (31) 

w r r 2 40 6 r 6 v ' 

S(r) = 0{r- w ), (32) 

^=-3^ + ^- < 33) 

These asymptotic solutions are later used as initial approximations for the numerical 
integration. 



3.4 Units 

In the system of units that we work G* = c = f 2 - = 1, where fio is the magnetic 
constant. In this system 

[P] = m; [M] = m; [m*] = m~ ; [6] = m~ 2 . 

In order to obtain dimensionless quantities we can use the fact that the differential 
equations system (j2"U|) - (l2"2"|) is invariant under the rigid re-scaling r — > Ar, m — > Am, 
P — ► AP and 6 — > A _2 6 where < A < oo. Therefore, given a solution to (|2"U|) - (j2"2"|) 
with one set of physical parameters (r h , M, P, m*,6, T), the rigid re-scaling produces 
new solutions with parameters (Ar^, AM, AP, X~ 1 m 3f , A _2 6, A _1 T). Here T denotes the 
temperature of the horizon. The dimensionless parameters can be obtained through 
such re-scaling if [A] = m~ x . Here, we introduce scale by choosing A = y/b. We will 
keep the notation for the dimensionless quantities unchanged. 
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4 Extremal solutions 



The massive-scalar-field black-hole solution has a much richer causal structure com- 
pared to the case with massless scalar field, since in the former case the presence of inner 
and extremal horizons is allowed. For a two-fold degenerate horizon from equations 
f j2TT) and (122jl having in mind that on the extremal horizon 



f(r e 







dm(r e 
dr 



1 
2 



hold, we obtain 



[V(ip e )-A\ip e )L(X t 



r 2 e - MVe)A\^ e ) [L(X e ) - X e d x L(X e )] > = 0. 



(34) 

(35) 
(36) 



Here and forth, unlike the previous section, r e will denote the extremal horizon. Again, 
subscript "e" denotes the value of the functions evaluated at r e . Equations fl35|) - fl36|) 
can be solved with respect to r\ 



-F x ± yjFl - 4a(y? e )F 2 



where 



Ft 
F 2 



dV{ip e 



46 



dip 
dip 



2F 2 



4:a(ip e )V(ip e ) 



dV(<f e 

dip 



(37) 

(38) 
(39) 



The root for which < r 2 e < oo for all values of ip e should be chosen. Then we 
substitute it in equation (1331) and obtain the following sophisticated non-linear equation 
for ip e 



1 _ -F 1 + ^F?-4aMF 2 ^ 



2F 2 

2V(ip e )-4bA 4 (ip e 



\ 



p2 

2b 



A\ip> e 



2F 9 



Fi + ^JFI - Aa(ip e )F 2 



(40) 



which we treat numerically for the STT with a(ip) = const = a. The solution for ip e 
is substituted back in (!37|) which, on its turn, gives us the radius of the degenerate 
horizons. Equation (140!) may have two roots for ip e which give two roots for r 2 e when 
substituted back in (!37|) . Which of them describes a degenerate event horizon? Using 
equation (fT5"|) we can find the value f" 



fl = -4 {V(<p e ) + A 4 (ip e ) [2X e d x L(X e ) - L(X e )]} 



(41) 
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Hence, we can use the fact that at the degenerate event horizon f" > to obtain a 
criterium which tells us whether the roots for (f e and r e obtained from (1401) and (157)) 
correspond to a degenerate event horizon. The roots for which f" < could correspond 
to an internal degenerate horizon resulting from the merger of two internal horizons if 
they exist at all. 

Equation (HOl) for ip e has two roots. The corresponding values r ei and r e2 are 
given in figure (CD). It should be noted that the analysis for existence of degenerate 
horizons does not take in consideration the boundary conditions. In other words, it is 
possible that for some values of the parameters the extremal black-hole solutions to be 
asymptotically non-flat. 

From the graphics we see that for fixed values of a and m*, a critical value of the 
magnetic charge P cr it exists, such that for P > P cri t solutions with degenerate event 
horizons exist and for P < P crit solutions with degenerate event horizons do not exist. 
For example, when a = 0.01 and ml = 0.8, P cr it ~ 0.38. A degenerate horizon is 
usually formed when for some values of the black-hole mass M two regular horizons 
merge. So, it is likely to expect that for P > P cr it black holes with two regular horizons 
exist for some values of M. On the contrary, when P < P crit we should not expect the 
existence of solution with more than one horizon. Our further numerical investigation 
confirms the expectations for the case P > P cr it- We should also note that for fixed 
P > P cr i t with the decrease of M an extremal solution is reached while for P < P crit 
— a naked singularity. 

The P — f" dependence is presented in figure fl2]). As we can see, for the whole 
interval of its existence the root r ei cannot be excluded as a candidate for a degenerate 
event horizon. The situation with the second root is more complicated. For most 
values of P, r e2 cannot correspond to a degenerate event horizon. For a small interval 
close to Pcrit, however, f" is positive. There are several ways to interpret r e2 in that 
interval. It could correspond to a degenerate internal horizon where the function / has 
a minimum. Since this analysis does not take into account the boundary conditions, 
another possibility would be that this root corresponds to a degenerate event horizon 
of a black hole with a different asymptotic behavior. The point where f'J turns to 
zero could correspond to a triply degenerate event horizon. Unfortunately we have not 
been able to check any of these speculations since the numerical simulations fail in that 
region of the parameter space. 

Alternatively, we could keep the values of a and P fixed and study the presence of 
degenerate event horizons for different m*. The numerical solutions for that case are 
shown in figure Q. Again, a critical value for scalar-field mass appears m* jCr it> such 
that for m* > m* jC rit solutions with degenerate event horizons exist and for m* < m* )C rit 
no solutions with degenerate event horizons exist. For example, when a = 0.01 and 
P = 0.4, m^crit ~ 0.26. Critical mass m* iC rit decreases with the increase of the magnetic 
charge P. 

The m* — f" dependence is presented in figure (@|. The comments about the two 
roots are the same as in the case presented in figure fl2J). 

The critical value P cr it depends on m* and a. In order to determine this dependence 
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Figure 1: Solutions for the positions of the degenerate horizons as a function of the 
magnetic charge P of the black hole. 
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Figure 2: Derivative f" evaluated on r e \ and r e 2 as a function of the magnetic charge 
P of the black hole. On a degenerate event horizon f" > must hold. 
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Figure 3: Solutions for the positions of the degenerate horizons as a function of the 
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Figure 4: Derivative f" evaluated on r e \ and r e 2 as a function the scalar field mass m*. 
On a degenerate event horizon f" > must hold. 



we will use again equations (l35|) and (I36l) . Let us introduce a new variable 



The variable k can be expressed as a function of <^ e from ( |36i) . Using (I3.1I) and ( 1421) we 
can express r 2 in the following way 

Then, in (|35|) we express I/(X e ) with k and from (j4"3"l) . Introducing a new variable 
y = «v9 e we obtain 



l = V2|P| ^q±gpi = V2|P|^„). (44) 



where 

2 



7 



4a 2 ' 

In f l44l) there are two free parameters \P\ and 7. Let us fix the value of 7. Then F(y, 7) 
has one maximum (y < 0). Depending on \P\ the RHS of (j44j) can be greater, equal 
to or less than the LHS and the equations have, respectively, two, one, or no solutions. 
So P cr itEI is the value for which (j44"j) has only one solution. The 7 — P cr i t critical curve 
is plotted in figure (jSJ). 



5 Numerical integration 

The nonlinear system (l2"01) - (122l) is inextricably coupled. Our aim is to obtain numerical 
solutions that describe asymptotically flat black holes. We split the problem in two 
boundary sub-problems - in the exterior and in the interior region of the black hole, 
which we solve subsequently. 

3 P is non-negative. 
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Figure 5: The 7 — P crit critical curve 

5.1 External problem 

First, the exterior region. For that region we can formulate a boundary value problem 
(BVP) for fl20])-([22]) with the boundary conditions 

lim mir) = M (M is the mass of the black hole in the Einstein frame), 

r— >oo 

lim S(r) = lim (fi(r) = 0. 
for the right-hand boundary at spatial infinity and 

f(r H ) = 

on the horizon. The left-hand boundary, namely the event horizon, is a priori unknown. 
Such BVPs are known in mathematical physics as BVPs of Stefan kind (see, for example 
|58j). For the location of the event horizon an additional condition is needed. It can 
come from the requirement all functions to be regular on the event horizon 

+ {4a(<p)A*(<p)[Xd x L(X) - L(X)]}\ r=rH . 

r=r H 

We can also think of the so formulated BVP as a non-linear analogue of a spectral 
problem with regard to parameter r#. 

Having in mind these features of the above posed BVP we treat it by using the Con- 
tinuous Analog of Newton Method (see, for example [19] , [50] , [13] ) . After an appropriate 
linearization the original BVP is rendered to solving standard vector two-point BVPs. 
On a discrete level almost diagonal linear algebraic systems with regard to increments 
of sought functions S(r), m(r), and (p(r) have to be inverted. 

5.2 Internal problem 

The exterior solutions can be continued inwards. The values of the functions and their 
derivatives on the event horizon are obtained a priori in the exterior problem. This 



df dip 
dr dr 



r=ru 



dV{(p) 
dip 
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allows an Initial- Value Problem (I VP) for the same system (l20l) - (!22l) to be formulated 
in the interior region r < r# a posteriori. The event horizon, however, is a singular 
point for the equation of the scalar field (1221) since the coefficient in front of ip" turns 
to zero there (/(?"#•) = 0) and the equation loses its leading order term. So, to pose a 
regular IVP we shift the initial point r# by small enough e > and choose for initial 
point th — £ instead r#. On the other hand the functions in question are smooth in 
the interval (r# — e, rn) and hence the following series expansions hold 



A similar shift is made also at the every inner horizon (if such is reached). The latter 
admits an algorithmic sequence of IVPs for finding possible inner horizons. For the 
numerical treating of the above posed IVP again the Continuous Analog of Newton 
Method is used. 

5.3 Some results 

Let us first consider the case P > P C rit- The numerical investigation of the solutions 
shows that for fixed values of the parameters a and m* the general structure of the 
solutions depends strongly on the charge-to- mass ratio P/M. For low enough P/M the 
obtained black-hole solutions have a single horizon, namely the event horizon. For high 
enough P/M (for P/M close to 1) the black holes have two horizons or one degenerate 
horizon (an extremal solution is reached when we decrease the mass M and keep P 
fixed). 

Two examples of solutions with a single horizon are shown (the functions /, 5 
and ip, respectively) in figures (jSJ), © and jS]). The values of the parameters in the 
presented examples are a = 0.01, ml = 0.8, P = 6.0 and two different masses of the 
black hole M = 21.5 and M = 25.0. As it can be seen in figure ([6]) for M = 21.5 the 
function / has negative extrema below the event horizon (the abscissa in that figure is 
in logarithmic scale). In figures ([7]) and (jSJ) the radius of the solution with M = 21.5 
is designated with "O" and for M = 25.0 - with "x." 

Solutions with two horizons and with one degenerate horizon are given in figures 
f fTUj) and flLTJ. The values of the parameters in these cases are a = 0.01, ml = 0.8, P = 
6.0 and black-hole masses M = 5.9804 for the extremal black hole and M = 8.0 for 
the black hole with two horizons. We can see that in accordance with our observations 
from Section [3] 5 (r) has an inflexion point in the extremum of ip(r). In figures (flQj) and 
(llip the radius of the extremal black hole is designated with "O" while the two radii 
of the solution with M = 8.0 are designated with "x"-es. 

The presence of numerical instabilities did not allow us to study thoroughly the 
interior region for the case P < P C nt- 



m(rn 
p(r H 




■- m(r H ) - m'(r H )e + o(e 2 ), 
S(r H )-6'(r H )e + o(e 2 ), 
^{rn) - p'{r H )e + o(e 2 ), 

$(r#) - <f>'(r H )e + o(e 2 ), where $(r) = (p'(r) 



(45) 
(46) 
(47) 
(48) 
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Figure 6: The metric function / for a = 0.01, ml = 0.8, P = 6.0 and two different 
masses of the black hole M = 21.5 and M = 25.0. The abscissa is in logarithmic scale. 



0,004- 



5{r) 



0,002- 



0,000- 



0,0 0,5 



M=21.5 
■M = 25.0 



-/A- 



1,0 40 

r 



-e- — j^- 



50 60 



Figure 7: The metric function 5 for a = 0.01, ml = 0.8, P 
masses of the black hole M = 21.5 and M = 25.0. 
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Figure 8: The scalar field ip for a = 0.01, = 0.8, P = 6.0 and two different masses 
of the black hole M = 21.5 and M = 25.0. 
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Figure 9: The metric function / for a = 0.01, ml = 0.8, P = 6.0 and two different 
masses of the black hole M = 5.9804 and M = 8.0. The abscissa is in logarithmic 
scale. 
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Figure 10: The metric function 5 for a = 0.01, ml = 0.8, P = 6.0 and two different 
masses of the black hole M = 5.9804 and M = 8.0. 
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Figure 11: The scalar field (p for a = 0.01, ml = 0.8, P = 6.0 and two different masses 
of the black hole M = 5.9804 and M = 8.0. 
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6 Thermodynamics 



Since the entropy of the black hole is related to the area, and respectively to the 
radius, of the event horizon in figures (fT2l) and (fTBl) we have given the M — vn and the 
P — th diagrams. In the presented cases a = 0.01, ml = 0.8. For these values of the 
parameters P cr it ~ 0.38. We have given two graphics - one for P = 6.0 > P cr it and one 
for P = 0.3 < P cr it- In the former case for masses in the interval M G [5.9804,21.1] 
the black holes have two regular horizons (in figure (1121) they are designated as event 
horizon and inner horizon) . The two horizons merge and an extremal black hole occurs 
at M = 5.9804. The value of event horizon obtained through numerical integration 
of the equations and the one obtained through solution of the algebraic problem in 
Section H] coincide within the approximation and round-off error, which is a very good 
test for the correctness of the obtained numerical results. In the latter case the black 
holes have a singe non-degenerate horizon. 

Let us try to obtain some intuition of the general behavior of our solutions through 
comparison to the solution of Reissner-Nordstrom (RN). In RN black holes exist for 
M > P. An extremal black hole occurs at P/M = 1. From the M — r# relation ( fT2l) 
we see that, similarly, in our case black-hole solutions exist for M > P. Unlike RN, 
however, in our case black holes exist also for M a bit less than P but still the extremal 
black hole occurs at P/M close to 1. So if we fix the value of M the maximal value 
of P is of the same order. We can define also M cr ; t as the mass of the extremal black 
hole when P = P cr i t . The relation P crit /M crit « 1 also holds. This dependence can be 
seen also in figure (fT3"|) where the P — r# relation is shown @. Let us set M = 0.2. 
Then according to our admission P cannot be much greater than 0.2 so it will be less 
than P cr it- Indeed, from figure (fTBl) we see that in that case a naked is reached when P 
approaches its maximal value (see the comments at the end of Section H]). For M = 3.0 
an extremal black hole appears when P approaches its maximal value. 




M 



Figure 12: The radii of the horizons of black holes as function of the mass M. For 
P = 6 the black hole has two horizon - event horizon and inner horizon. For P = 0.3 
the black holes have a single horizon - the event horizon. 



4 On that figure only the radii of the event horizons are shown. 
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Figure 13: The r# — P relation. 



A knowledge about the general stability of the obtained solutions can be obtained 
through the application of the so-called "turning point" method (we refer the reader 
to [5H [52J E3 El] for a detailed discussion on the "turning point" method and also to 
[55l H2J [3] for the application of the method to study the thermodynamical stability of 
black holes). According to that method, in micro-canonical ensemble 0, a change of the 
stability reveals itself on the M — T _1 diagram as bifurcation or turning points. Here 
"bifurcation point" is used to denote a point where brunching of equilibrium sequences 
occurs while a "turning point" is such point where two equilibrium sequences merge 
with a vertical tangent. The absence of such points on the M— T _1 diagram means that 
if at least one point on the equilibrium sequence is stable then the whole equilibrium 
sequence is stable. 

The M — T~ l diagram for the studied solutions is given in figure ( 1T41) . Again we 
present the two cases P = 6.0 > P cr ; t and P = 0.3 < P cr it- In both of them, when 
M is large T™ 1 « M, i.e., the solutions approach the Schwarzschild black hole. The 
stability of the Schwarzschild black hole within the theory of Brans-Dicke has already 
been proved [57] so we expect that for large M our solutions are stable. Since in both 
of the presented cases no turning or bifurcation point appears on the diagram we can 
expect that the entire equilibrium brunches are stable. With the decrease of M the 
inverse temperature T _1 of the solutions with P > P CTlt increases unboundedly since 
an extremal solution is reached, while for P < P CT it, T^ 1 goes to zero since an object 
with zero radius of the horizon, i.e. a naked singularity, is reached. In figure ffl5l) the 
P — T relation is shown. For M > M crit an extremal black hole is reached with the 
increase of P so the temperature of the horizon goes to zero. For M < M crit a naked 
singularity is reached with the increase of P and the temperature of the horizon rises 
unboundedly. 

5 We keep the magnetic charge P of the black hole fixed. 
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Figure 14: The M — T 1 relation. 
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Figure 15: The P — T relation. 
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7 Stability of the solutions against spherically sym- 
metric perturbations 

7.1 Stability of the magnetically charged black holes 

In order to examine the stability of the solutions against radial perturbations we apply 
the scheme presented in details in [3]. First, we establish the stability against ra- 
dial perturbations of the magnetically charged solutions. Then, study the electrically 
charged case through the electric-magnetic duality rotations. We use the following 
ansatz for the metric of a spherically symmetric, time dependent space-time 



ds 2 = -e^dt 2 + e x dr 2 + e\dd 2 + sin 2 



(49) 



where 7, \ and /3 are functions of r and t. Then, the perturbed the fields are presented 
in the following way 



7(r,t) 



-28{r) + In /(r)J + A 7 (r, t), x {r, t) = - In /(r) + A X {r, t), 

P(r,t)=2hir + Ap(r,t), <p(r,t) = <p(r) + A<p(r,t), (50) 



where 6(r), f(r) and ip(r) give the static background solution, and Aj(r,t), Ax(r,t), 
Af3(r,t) and Aip(r,t) are small time-dependent perturbations. Let us impose the 
convenient gauge A(3(r,t) = 0. Then simply reduces to r 2 . 

In the case of spherically symmetric perturbations that we consider the equations 
decouple and the system reduces to a single equation for the perturbations of the scalar 
field 

V (0) V (0)M A ^ _ u^ Aip = 0j ( 51 ) 

where 



U{r) 



2\ 1 - 2r' 



V-A\y)L{X) \ d r if(r 



n 2 



+ 4r 



d r <f{r) 



dV (fV_ 

dip dip 2 



+ A\ V ) l%rd r Lp(r)a{ip) + 4d 9 a((p) Xd x L{X) - L(X) 



- lQa 2 (ip)A\<p) X 2 d 2 x L(X) - Xd x L(X) + L(X) , (52) 

and is the co-derivative operator with respect to the static background. Since the 
background solution is static equation ( 15T1) admits a separation of the variables. Using 
the following substitution 

Aip(r,t) =ip(r)e iuJt , (53) 
from (IBTjl we obtain an equation for the spacial part of the perturbations 



/(r)e -*,- 



/(r)e -*,£ 



+ lf(r) e - 2 ^) d ± + ^ = f( r )e- 2S ^U(r)^, (54) 
r dr 



where uj 2 acts as a spectral parameter. Equation ([54]) can be cast in the form of 
Schrodinger equation. In order to do this we will use the tortoise radial coordinate 
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r* and will make a proper substitution for i/j 

dr 



dr* 



f(r)e-s<r) 

The following Schrodinger-like equation 

d 2 u(r*) 



i/j(r) 



u(r) 



dr 2 



+ uj u(r*) = U eS (r*)u(n) 



(55) 



(56) 



is reached. The effective potential is 

U cS (n) = f(n)e- 25 ^{u(n) - 2V(<p) + 2A\ V )L(X) + -±-[1 - /(r*)]}. (57) 

For r varying in the interval [r#, oo), where r# is the radius of the event horizon, the 
tortoise radial coordinate r* 6 (—00,00). From this point we can use the techniques 
from the standard quantum mechanics to study the properties of the small perturba- 
tions. 

In the numerically studied region of the parameter space the effective potential U e g 
is non-negative which means that the studied black holes are stable against spherically 
symmetric perturbations. Two illustrative cases are given in Figure (1161) for a = 0.01, 
ml = 0.8, P = 6.0 and two values of the black-hole mass: M = 10.0, for which the 
black hole has a regular event horizon, and M = 5.98 corresponding to an extremal 
black hole. Here we should note that in the figure the effective potential is presented 
in terms of the radial coordinate r. In terms of the tortoise coordinate r*, U e g would 
be more stretched but still non-negative. 
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Figure 16: The effective potential U e Q for a = 0.01, ml = 0.8, P = 6.0 and two values 
of the black-hole mass M = 10.0, for which the black hole has a regular event horizon, 
and for M = 5.98 corresponding to an extremal black hole. 



7.2 Stability of the electrically charged black holes 

The perturbations of the fields in the electrically charged case can be obtained from 
the magnetically charged case through the electric-magnetic duality rotations. These 
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transformations preserve the functions of the metric and the scalar field. The pertur- 
bations of the electromagnetic field, however, will be non-vanishing in the electrically 
charged case. For them the duality rotations give (see [5]) 

AI^=i(*GMA x + A7) 

+ Aa{^)Xd 2 x L{X){*F)^ A V - d x L(X)(*AF)^. (58) 

Thus, the perturbations of the electromagnetic field can be expressed in terms of the 
functions of the background, static solution and the perturbations of the magnetically 
charged solution. So the perturbations of the electrically charged solution will remain 
bounded with time as long as the perturbations of the magnetically charged solution 
are. 

8 Conclusion 

In the present work numerical solutions describing charged black holes coupled to 
non-linear electrodynamics in the scalar-tensor theories with massless scalar field were 
found. Since an electric-magnetic duality is present in the used electrodynamics, only 
purely magnetically case was studied here. For the Lagrangian of the non-linear electro- 
dynamics the truncated Born-Infeld Lagrangian was chosen and scalar-tensor theories 
with massive scalar field and positive coupling parameter were considered. As a result 
of the numerical and analytical investigations, some general properties of the solutions 
were found. The theory we considered admits the existence of extremal black-hole 
solutions unlike the case with massless scalar field. 

Acknowledgments 

This work was partially supported by the Bulgarian National Science Fund under 
Grants VUF-201/06, DO-02-257, Sofia University Research Fund No 074/2009 and 
Technical University Sofia Grant 091nil43-l 1/2009. 

References 

[1] I. Stefanov, S. Yazadjiev and M. Todorov, Phys. Rev. D 75, 084036 (2007). 

[2] I. Stefanov, S. Yazadjiev and M. Todorov, Mod. Phys. Lett. A 23 (34), 2915 

(2008) , arXiv: 0708.4141. 

[3] I. Stefanov, S. Yazadjiev and M. Todorov, Class. Quant. Grav. 26, 015006 

(2009) . 

[4] M. Born and L. Meld, Proc. R. Soc. London A143, 410 (1934). 
[5] E. Fradkin and A. Tseytlin, Phys. Lett. B 163, 123 (1985). 



23 



[6] E. Bergshoeff, E. Sezgin, C. Pope, and P. Townsend, Phys. Lett. B 188, 70 
(1987). 



[7 
[8 
[9 

[10 
[11 
[12 
[13 
[14 
[15 
[16 
[17 

[18 
[19 
[20 
[21 
[22 
[23 
[24 
[25 
[26 

[27 
[28 
[29 



R. Metsaev, M. Rahmanov, and A. Tseytlin, Phys. Lett. B 193, 207 (1987). 
R. Leigh, Mod. Phys. Lett. A4, 2767 (1989). 

A. Garcia, H. Salazar and J. Plebanski, Nuovo Cimento Soc. Ital. Fis. B 84, 65 
(1984). 

M. Demianski, Found. Phys. 16, 187 (1986). 

D. Wiltshire, Phys. Rev. D 38, 2445 (1988). 

H.P. Oliveira, Class. Quantum Grav. 11, 1469 (1994). 

G. Gibbons and D. Rasheed, Nucl. Phys. B 454, 185 (1995). 

G. Gibbons and D. Rasheed, Phys. Lett. B 476, 515 (1996). 

E. Ayon-Beato, A. Garcia, Phys. Rev. Lett. 80, 5056 (1998) 
T. Tamaki, T. Torii, Phys. Rev. D 62, 061501R (2000) 

M. Novello, S. Perez Bergliaffa and J. Salim, Class. Quantum Grav. 17, 3821 
(2000). 

K. Bronnikov, Phys. Rev. D 63, 044005 (2001). 

G. Gibbons and C. Herdeiro, Class. Quantum Grav. 18, 1677 (2001). 

H. Yajima and T. Tamaki, Phys. Rev. D 63, 064007 (2001). 

M. Gurses and O. Sarioglu, Class. Quantum Grav. 20, 351 (2003). 

A. Burinskii and S.R. Hildebrandt, Phys. Rev. D 65, 104017 (2002). 

S. Fernando and D. Krug, Gen. Relat. Grav. 35, 129 (2002). 

T. Dey, Phys. Lett. B595, 484 (2004). 

T. Kumar Dey, Phys. Lett. B 595, 484 (2004). 

I. Stefanov, S. Yazadjiev and M. Todorov, Mod. Phys. Lett. 22 (17), 1217 
(2007). 



D. Rasheed, arXiv: hep-th/9702087 



N. Breton, arXiv: gr-qc/0109022, 



N. Breton, arXiv: hep-th/0702008 



24 



[30] R. Cai, D. Pang and A. Wang, Phys. Rev. D70, 124034 (2004), arXiv: 



hep-th/0410158 



[31 
[32 
[33 
[34 
[35 
[36 
[37; 

[38 
[39 
[40 
[41 
[42 
[43 

[44; 
[45 
[46 
[47; 
[48 

[49 
[50 

[51 



T. Tamaki, J CAP 05, 004 (2004). 

G. Clement and D. Gal'tsov, Phys. Rev. D 62, 124013 (2000). 
T. Tamaki, T. Torii, Phys. Rev. D 64, 024027 (2001). 
S. Yazadjiev, Phys. Rev. D 72, 044006 (2005). 

A. Sheykhi, N. Riazi and M. H. Mahzoon, Phys. Rev. D 74, 044025 (2006). 
A. Sheykhi, N. Riazi, Phys. Rev. D 75, 024021 (2007). 

T W. Chemissany, Mees de Roo, S. Panda, Class. Quant. Grav. 25, 225009 
(2008). 

A. Sheykhi, Phys. Lett. B 662, 7 (2008). 
A. Sheykhi, IJMPD (2008) to appear, arXiv: 0801.4112. 
R. Gregory, J. A. Harvey, Phys. Rev. D 47, 2411 (1993). 
J. Home, G. Horowitz, Nucl.Phys. B 399, 169 (1993). 
T. Tamaki,Phys.Rev. D 66, 104021 (2002). 

S. Yazadjiev, P. Fiziev, T. Boyadjiev and M. Todorov, Mod. Phys. Lett. A 
16, 2143 (2001). 

N. Banerjee and S. Sen, Phys. Rev. D47, 104024 (1998). 

Th. Sotiriou, V. Faraoni, St. Liberati, Int.J.Mod.Phys. D 17, 399 (2008). 

V. Faraoni, Sh. Nadeau, Phys. Rev. D 75, 023501 (2007). 

Th. Helbig, Ap.J. 382, 223 (1991). 

C. M. Will, The Confrontation between General Relativity and Exper- 
iment, Living Rev. Relativity 9, (2006), URL (cited on 29.01.2008): 
|http: //www.livingreviews.org/lrr-2006-3[ |arXiv:gr-qc/0103036| 

M. K. Gavurin, Izvestia VUZ, Matematika 14(6), 18-31 (1958) (in Russian) (see 
also Math.Rev. 25(2), #1380 (1963)). 

E. P. Jidkov, G. I. Makarenko and I. V. Puzynin, in Physics of Elementary 
Particles and Atomic Nuclei (JINR, Dubna, 1973) vol. 4, part I, 127-166 (in 
Russian), English translation: American Institute of Physics, p. 53. 

J. Katz, Mon. Not. R. Astr. Soc. 183, 765 (1978). 



25 



[52] R. Sorkin, Astrophys. J. 249, 254 (1981). 

[53] R. D. Sorkin, Astrophys. J. 257, 847 (1982). 

[54] G. Arcioni, E. Lozano-Tellechea, Phys. Rev. D72, 104021 (2005). 

[55] J. Katz, Mon. Not. R. Astr. Soc. 189, 817 (1979). 

[56] T. Torii, T. Tamaki, K. Maeda, Phys. Rev. D68, 024028 (2003). 

[57] O. J. Kwon, Y. D. Kim, Y. S. Myung, B. H. Cho, Y. J. Park, 
Phys. Rev. D34, 333 (1986). 

[58] A. N. Tikhonov, A. A. Samarskii, Equations of Mathematical Physics, 1990, 
Nauka, Moscow, (in Russian) (English translation: Dover Publications, NY, 
1990). 



26 



